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Abstract 

A detailed analysis of the remainder obtained by truncating the Euler series up to the 
nth-order term is presented. In particular, by using an approach recently proposed 
by Weniger, asymptotic expansions of the remainder, both in inverse powers and in 
inverse rising factorials of n, are obtained. It is found that the corresponding expanding 
coefficients are expressed, in closed form, in terms of exponential polynomials, well 
known in combinatorics, and in terms of associated Laguerre polynomials, respectively. 
A study of the divergence and/or of the convergence of the above expansions is also 
carried out for positive values of the Euler series argument. 

1 Introduction 

We consider the so-called Euler series (ES henceforth), defined as 

oo 
m— 

where z is a nonnegative number (possibly complex). The series in Eq. ([TJ has 
a zero convergence radius and gives a coded representation of the function 

£(z) = e -^^E 1 ( 1 -), (2) 



z \z / 

where E\{x) denotes the exponential integral, denned by pQ 



e 1 (x)= I ^pLi) dt; (3) 



for | argai| < it. The function in Eq. ([2]) is called the antilimit of the series in 
Eq. ([T]) and, according to Euler, can be viewed as the mathematical entity whose 
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expansion gives rise to the divergent seriesQ Incidentally, the ES represents the 
paradigm for many other factorially divergent asymptotic inverse power series 
occurring in special function theory, or arising from large-order perturbation 
expansions in theoretical physics. 

The resummation of factorially diverging series like the ES can be success- 
fully operated through the use of different strategies [TU]. Among them, those 
based on nonlinear sequence transformations have proved, especially in recent 
times, to achieve the retrieving action in an effective way. The common feature 
of several types of sequence transformations is the following decomposition of 
the nth-order partial sum s n — X)L-=o a k °f the starting series: 

s n = s + r n , (4) 

where s denotes the antilimit and r n = — Xm+i a k the nth-order remainder. 
For the ES the decomposition in Eq. ((4]) can be straightforwardly derived by 
writing 

n 

Sn = (-z) m m\, (5) 

m— 

and by expanding the factorial as m! = / °° dt exp(— t)t m , thus obtaining 
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dt e -^tH-(- z y^ fdt exp (-^" +1 . 

l + Zt V ' Jn l + Zt 



where use has been made of the explicit expression of the partial sum of the 
geometric series, i.e., 

1 - x n+1 

' ' ■ (7) 



E 



1 ~ x 

m—0 

On comparing Eq. ([6]) with Eqs. ((2J and (J4j), it is seen at once that 
s = £(z), 

ex p(-t) t n+1 (8) 
1 + zt 



r n = -(-z) n+i / dt 



The development of nonlinear sequence transformations is connected to the 
theory of converging factors [TUl According to it, the remainder r„ is 
expressed as the product between the first term of the series not included in the 



There is a nice quotation of a letter from Euler to Goldbach (1745) in Ref. |10l . which 
reads: "Summa cuiusque seriei est valor expressionis illius finitae, ex cuius evolutione ilia 
series oritur, " which may be translated from the Latin language as: "The sum of any given 
series is the value of the specific finite expression whose expansion gave rise to that same 
series." 
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partial sum, a n +i, and a converging factor, say ip n , chosen in such a way that 
the relationship 

s n — s + a n+ iip n: (9) 

is satisfied. In particular, from Eq. © it follows that, for the ES, an integral 
representation of the nth-order converging factor (or terminant) is (see Ref. [T2] , 
Ch. 21) 

1 f 00 , exp(-t)t n+1 
(n + 1)! J l + zt 

The aim of the present work is to find expansions, both in inverse powers 
and in inverse rising factorials of n, of the converging factor of the ES in Eq. (JTJ) 
in such a way that the nth-order remainder in Eq. flSJ can be expressed through 
the forms 



a 



n+i 7 — ; rr, (H) 



fc=o 



and 



oo 

\ " 



r n = a n+1 > - — ■ — — , (12) 

where a denotes a positive parameter, (•)& denotes the Pochhammer symbol, 
defined by 

r(a + fc) 

Wk = r , s , (13) 

r(a) 

and where the infinite sequences {ck} and {dk} are independent of n. In par- 
ticular, the last prescription revealed crucial for the Levin |15] and the 5 (or 
Weniger) [50] transformations, which are strictly related to Eqs. (fTTjl and (fT2"]) , 
respectively, to be derivedH On the other hand, it should also be noted that 
different types of asymptotic expansions of the converging factors for the expo- 
nential integral function have already been found in the past \17\ [5] , for instance 
by using the Airey's approach [2]. 

To find the closed-form expressions of the sequences {cfc} and {dk} we shall 
use the approach recently introduced by Weniger in Ref. [35] • According to it, 
the starting point is the decomposition scheme in Eq. ((4]) , from which it follows 
that the truncation error r n associated to the partial sum s n must satisfy the 
first-order difference equation 

Ar„ = a n+1 , (14) 

or 

r " +1 = r " = 1, (15) 

where A denotes the forward difference operator with respect to n, i.e., such 
that A/(n) = f(n + 1) — /(n). In Ref. [22] such approach was used to reproduce 



2 For an extensive review about applications of the Weniger transformation, see for instance 
Refs. 10 21 6 7 and references therein. 
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the Euler-Maclaurin formula for the remainder of the Dirichlet series for the 
Riemann zeta function. In the case of the ES, the same approach led to the 
exact expression of the first few terms of expansions similar to those in Eqs. 
and (fT2")l . In the present work we will show that the whole sequences {ck} and 
{dk} appearing in Eqs. (|TT|) and (fT2")l . respectively, can be obtained, for the case 
a = 1, through simple, analytical closed-form expressions. In particular, what 
we will find is that the c^'s are expressed by exponential polynomials [4l [9l HTj. 
while the coefficients dk turn out to be proportional to associated Laguerre 
polynomials pQ. 



2 Theoretical Analysis 

2.1 The Weniger approach for building up asymptotic 
expansions of truncation errors 

As anticipated in the previous section, the approach proposed by Weniger in 
Ref. [22] will be pursued to express, in closed-form terms, the two sequences 
{cfc} and {dk} in Eqs. (fTTj) and (fl"2"|) . for a — 1. We begin with the asymptotic 
inverse power series, needed for the Levin transformation. The inverse factorial 
expansion, related to the Weniger transformation, will be subsequently derived 
starting from the former. However, to obtain the asymptotic expansion of the 
converging factor as in Eq. (TTT|) . following the prescriptions given in Ref. [2"2"] . 
an intermediate step is necessary. It consists in replacing, in Eq. (|15[) . the 
remainder r n by the estimate, say r n m \ given by 

m m 

rn - °» £ j^Tjk = (-*) E (^TIp (16) 

fc=0 v ; fc=0 v ; 

where {70, 71, ... , 7 m } are unknowns quantities which must be indipendent of n. 
It should be noted that in Eq. (fT6"|) the truncation error of the ES is represented 
as the last term included in the partial sum multiplied by a truncated inverse 
power series in n + 1. The correct converging factor will be then obtained 
starting from the knowledge of the 7fc's. The key point in the approach of 
Ref. [2 2) consists in substituting from Eq. (|16[) into Eq. (II 5[) and in requiring 
that the subsequent equation be satisfied up to the power n _m , i.e., that 

(m) _ (m) 

" +1 r " =l + 0{n- m - 1 ), (17) 

O-n+l 

for n — > 00. After some algebra, it is obtained (2"2"| 



(m) (m) m m 

C+; - a; 7fe 7fc 



(-l/a;)"+ :L (n + 1)! n+l^(n + l) fc ^ (n + 2)' 



= 1 + 0{n- m - x ), 
(18) 
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where, for convenience, it has been set x = 1/z. Note that Eq. (fT5|) is asymp- 
totically equivalent to 

mm m ? 

x ST Jk 7fc «fe,0 -m-l\ / 1n s 

k=0 k=0 K ' k=0 

which led to a linear system for the coefficients jk that was explicitely solved 
by Weniger in Ref. [22] for m = 4. 

We are now going to prove that such system can be solved in closed form 
for any values of to. To show this, we start by writing 

tci tci m+1 m 

fe=0 fc=o v y fe=l fe=l v ; 

and we note that 



i i i _ 1 v ^ 

7k 



13 



(n + l) fe n fe (l + l/n) fc n fc nJ j\ 

= f t^L ( fc + i ~ 9' = f (=1)^ + = (21) 

1 ^ ^ T*3 (fc-1)! ^ ! Vfc-1/ 

Furthermore, on substituting from Eq. (I2T1) into Eq. ([20)1, we have 

m m 

n 2^ n k 2.^ ( n+ l)k 

k=0 fc=0 v ' 



fc=l fc=l j=fc 



fc— 1 fc— 1 j—k 



x 7m \ - r _ n fc+i Tfef^ 1 

n m+i T- V J U - 1 

fe=l j=m+l v 



m mm / ■ -i \ 

fe=i fe=i j=fe v 7 
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or, by interchanging the symbols j and k, 



-E^ + E 



7fc 



n * — ' n 

k=0 



k=0 



to+E ^ 



fc=i 



(n + 1) A 



fc - 1 
j-1 



7j 



(23) 



•0(n 



-m—l\ 



Equation (f2"3"l) . together with Eq. (|19p . leads to the following linear system for 
the 7fe's: 

7o = 1, 



k 

3 = 1 



k-1 



3 



(24) 



, 7i=0. 



The solution of such system can be expressed in closed form simply by evaluating 
the quantity 

f n — 1 N 



-x 



E 

k=l 



k-1 



7fe-i> 



(25) 



which, by taking Eq. (|24|) into account, takes on the form 

71 fe 



E 

fe=i 



n — 1 

fc - 1 



7.- 1 = EE(- 1 ) fe+j 

fc=l j=l 



n - 1\ (k-1 
k-l)\j-l 



1i 



(26) 



j=l k=j 



n-l\ fk-1 



k-1) \3 



! =E(-!) 2J '^' 



3 = 1 



where formula 4.2.4.45 of [18] has been used. It is worth exploring Eq. (|26|) in a 
deeper way. Actually, this equation already contains the closed-form expression 
of the 7„ coefficients, as we shall see in a moment. 



2.2 The asymptotic expansion of the truncation error and 
the exponential polynomials 

Due to its importance, we rewrite Eq. (f2l))) as 




It must be noted that there exist a whole class of functions satisfying the relation 
in Eq. Q27p. Such functions are called exponential, or Bell, polynomials [J], and 
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will be denoted <^ n (a;)o These polynomials are denned through the following 
generating function formula: 

00 rn 

<Mz)^=exp{z[exp(i)-l]}, (28) 

n=a 

and have the explicit expansion 

n 

<t )n {x)=Y J S{n,k)x\ (29) 

fc=0 

where S(n, k) denotes the Stirling number of the second kind, which is defined 
throughj 

n 

Z " = (-l)«^(-l)*%,fc)( 2 ), (30) 

fc=0 

Exponential polynomials satisfy the recurrence relationship 

4>n+l(x) = x[(j>' n (x) + <p n (x)], (31) 

with 4>q(x) — 1 and, more importantly, they fulfill the relation 

<t> n (x) =xJ2 (lZ]) <f>k-i(x), (32) 



fe=i 



which, when compared to Eq. (f2~Tl) , shows that the coefficient j). is proportional 
to cj>k(—x) and, by virtue of the initial condition 70 ~ 1, that 

lk=4>k{-x). (33) 

Equation (|33j) represents one of the main results of the present work. According 
to it, we shall express the nth-order remainder of the ES through the following 
asymptotic series: 

which, as pointed out at the beginning of Sec. 12.11 is not yet of the desired form 
given in Eq. (fTTj) . To find the correct sequence {ck} for a = 1, it is sufficient to 
recast Eq. (|34|) as 

_ Qn+l ^ 7fc _ -(/) k - 1 (-l/z)/z 

rn ~ z(n+l) ^ {n + l) k £ {n + lf ' [6b) 



3 Here and in the following we are going to use, for the Bell polynomials, the notation given 
in the recent review by Boyadzhiev [9]. Furthermore, note that the numerical evaluation of 
Bell polynomials 4> n {x) of arbitrary order is currently implemented, up to arbitrary precision, 
whithin the Mathematica platform through the command BellB[n,x]. 

4 Here and in the rest of the paper we use, for the Stirling numbers of the first and of the 
second kind, the notation adopted in Ref. 
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which, once compared to Eq. (jlip . leads to 

0, k = 0, 

Ck = { a, ( i / ^ ( 36 ) 
-^1™, k>0. 

z 

2.3 The factorial expansion 

The transformation of an inverse power series to a factorial series can be ac- 
complished with the help of the Stirling numbers of the first kind, say s(n, j), 
which are defined through 

n 

(-l) n (z) n = J2^,j)(-z) j . (37) 

In particular, the Stirling numbers of the first kind occur in the factorial series 
expansion of an inverse power, namely |22] 

JL = (_!)* y fCj = y- ( _iy gCj 

C fc h (Oj ft (Oj ' 

(38) 

valid for fc > 1, where in the last passage use has been made of the fact that 
s(n,m) = when m > n. Accordingly, as pointed out in Ref. |22) . given an 
asymptotic power series of the form Y2k=o C k/Q ■> the following identity can be 
established: 

oo OO , 



where 



c , fc = 0, 



d k = { k (40) 

(-^^(-i^^-i.i-iJci, *>1- 

Then, on substituting from Eq. (13T>1) into Eq. (T^U|) , after some algebra it is found 
that the expanding coefficients {dk} in Eq. (|12[) are given by 



z 

where the function i()k(x) is defined as 

k 




= < , (41) 

fc > 0, 



^(i) = (-l) fc ]T (-1)' 8{k,j)^(x). (42) 



3 The ip n (x) polynomials 



On substituting from Eq. (|4"Tj) into Eq. (|12p . we eventually obtain 

-^fc-i(-l/2)A 



^ -y fc _ U -i/^/^ 



In the next section it will be proved that tp n (x) is proportional to the associated 
Laguerre polynomial of orders n and -1 [TJ. 

3 The ip n (x) polynomials 

From the definition given in Eq. (1421) . on using Eq. (|29[) we have 

n j 
^ n {x) = (-1)" J2(-iy S (n,j) S(j,k)x k = 

j=0 k=0 
n n 

= (-1)" xk E s (^') = ( 44 ) 

n n 

= (-l) n ^2x k ^2(-iys(n,j)S(j,k), 
k=0 3=0 

where in the last passage use has been made of the fact that S(j, k) = for 
k > j. The expanding coefficients in Eq. (|44|) can be given a closed form. To 
show this, we first recall formula 24.I.4.C of Ref. [T], i.e., 



^ = ^£(-W*W (45) 



which gives 



£ (-1)'' s(n, j) S(j, fc) = £ (-1/ m £ S (n, j) = 

(46) 

(— l) fe +'" 



EH)' ; Wn, 



fc! 

£=0 

where in the last passage the definition of the generating function of the Stirling 
numbers of the first kind, given in Eq. ()37|) . has been used. Furthermore, on 
taking into account that 
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Eq. (|4"6"]l becomes 

±(-iy S (n, ] )S(j,k) = (-ir^±(-iy( k ) + (48) 
i=o " £=o \ / \ / 

and, by using formula 4.2.5.26 of Ref. [18], it is obtained 

On substituting from Eq. (|4"9[) into Eq. ([44)1. the polynomial expansion of ip n (x) 
turns out to be, for n > 0, 

n k / \ 

^) = (-i)'E^ I ) I Qj ! (so) 

from which it follows that 

i> n (x)=n\L< h -V(-x), (51) 

where L^(-) denotes the associated Laguerre polynomial of orders n and a pQ. 
Finally, on substituting from Eq. (|51l) into Eq. (|4"3")l . the factorial expansion of 
the nth-order remainder of the ES reads 



kl 

-^E^^W^. (52) 



fc! 

which, together with Eq. (|3"5)) , constitutes the main result of the present paper. 

4 Discussions 
4.1 Preliminaries 

First of all, we note that the remainder r n in Eq. (j8]) can be evaluated in closed 
form, and turns out to be given by 

6XP (z) /1\ 
r n = -a„+i E n+2 {-)> (53) 

which, once compared with Eqs. (|3"5j) and ([521 gives at once the following two 
expansions for the exponential integral function: 



^«+i(ar) = 2^ , (54) 

fc=0 
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and 

E n +i{x) = V -—— L^ix) exp(-x), (55) 
feu (n)fe+1 

respectively. Moreover, Eqs. (|54|) and (1551 can be cast, by using the connection 
between the exponential integral and the incomplete gamma functions pQ, in 
the following form: 

r(-n, a;) = ar" exp(-a;) (f> k (-x), (56) 

and 

T(-n, x) = x~ n exp(-o.) V — — (x), (57) 

respectively. In particular, the expansion in Eq. (151)1, displays a structure very 
similar to an asymptotic expansion recently found in Ref. [5] for the (lower) 
incomplete gamma function 7(A, x) [T], with A > 0. 



4.2 Analysis of the convergence of the two series for positive 



x 



It is worth studying the character of the series in Eqs. (|5"4")) and (|55|) . for x > 0. 
As far as the series in Eq. (l54l) is concerned, we use the asymptotics, for k ^> 1, 
of the Bell polynomials, recently reviewed in Ref. [13] ■ In particular, for x > 
we have [13] 



exp(-x) 0fe(-x) ps fc!^/ — 



exp < fc 



sin w sm <z> 
log cos w 

If If 



sm ip 



COS </? 



sin 2 if 



1/4 



(58) 



x sm 



k ( 7T — If + 



sin 2 



where 



7T 

2 



arccos 



1 — if cot y 



1 V4 



(1 — f cot (p) 2 + sin if 



and where tp is the solution of the equation 



(59) 



x 



sm if 



exp 



9 



sin ip 



cos y 



(60) 



For any fixed value of x, we are interested in estimating the behavior of exp(— x)<f>k{— x) 
for k — y oo. In such limit the solution of Eq. (|6"0"|) tends to 7r~. In the same 
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limit, rj(ip) it/2, so that the whole term sin 



k I 7T — ip + 



be replaced by 1. Secondly, in the limit ip — > ir , we have 



simp 



9 



sin if 



cos tp 



sin 2 <p 



1/4 



sin 2 ip 



1/2 



and 



can 



(61) 



exp < fc 



sin G3 sin <p 
log cos ip 



sin </? 



exp fc 



sin 99 



so that 



exp(— x) 4>k{—x) ~ kl 



sin </? 



fe+1/2 



exp fc 



sin 99 



9 



(62) 



(63) 



which, once inserted into Eq. (|54j) . shows that the asymptotic series displays a 
factorial divergence. 

We now prove that the factorial series in the r.h.s. of Eq. ([55]) is, for x > 0, 
convergent. This can be done by noting that, for n > 1, 



fc! 



1 (fc + 1)! 



1 



1 x 2 x . . . (fc + 1) 



< 



1 



(n) fe+ i fc + 1 (n) fe+ i fc+1 nx (n+1) x ...(n + fc) - fc+l : 

(-1), 



(64) 



which, by taking into acount the asymptotics of L k (x) for large fc, namely, 



L[7 1 \x)exp(~x)K — cxp(-x/2)x 1 ^k~ 3 ^ cos (2^ + J 



(65) 



leads, for sufficiently high values of fc, to 



fc! U~ k X \x) exp(— x) 



(n)k+i 



< 



H 1] (x) exp(-x) 



(66) 



< 1 exp(-x/2) z 1/4 1 ^, exp(-a:/2) x 1 / 4 1 



fc + 1 v^F fc 3 / 4 fc 7 /4 : 
that proves the absolute convergence of the factorial series. 



4.3 Some remarks for negative values of the ES argument 

If z < 0, the ES in Eq. ([TJ) becomes a nonalternating, divergent asymptotic power 
series in z, which has been used in the literature as a paradigmatic example of a 
series that cannot be resummed by the Levin and Weniger transformations [14] . 
The reason for such inability in the resummation process is strictly related to the 
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fact that the negative real axis coincides with the branch cut of the exponential 
integral function E\(l/z). As a matter of fact, the (all positive) single terms of 
the ES in Eq. ([1} cannot reproduce the imaginary part of the function E\(l/z) 
which, for z < 0, is evaluated as [T] 

E x (-±i0) =-Ei(-i ) TiTr, (67) 



, z / \ z 

where Ei(-) denotes the function pQ 



Ei(x)=V I GXP( ^ dt, x>0, (68) 



t 

with the symbol V I ... denoting the principal value operator, i.e., 



,--J&U_- + /. ■■■)■ m 

It must be stressed that, for z < 0, also the asymptotic and the factorial series in 
Eqs. (|35|) and ()52j) become nonalternating, their single terms being all positive. 
For the former this can be proved by noting that the exponential polynomials 
4>k are strictly positive, for any k, when their argument is positive, [this can be 
seen from the recurrence relation in Eq. (|32[)]. For the latter, the nonalternating 
character directly follows from the expansion in Eq. (|50|) . On the other hand, 
from the integral representation in Eq. (jiop . it appears that also the converging 
factor tp n of the ES presents a branch cut coincident with the negative real axis 
z < 0. Accordingly, we conclude by saying that for those cases in which the ES 
cannot be successfully resummed by Levin and Weniger transformations, also 
the series expansions of r„ in Eqs. (pl5j) and (|52"j) seem to lack their effectiveness 
in giving a meaningful representation of the remainder. 



5 Conclusions 



The understanding of the retrieving action, as well as the development of new 
types of nonlinear sequence transformations aimed at resumming different classed 
of divergent series requires the large index asymptotics of the corresponding 
truncation errors to be investigated. A general approach for achieving such task 
has recently been proposed in Ref. [22]. In the present paper we showed that, 
for the (factorially divergent) Euler series, such approach allows the nth-order 
remainder to be represented via asymptotic and factorial expansions involving 
exponential and associated Laguerre polynomials, respectively. The convergence 
of the above expansions has also been investigated. 
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